Atomistic mechanisms and diameter selection during nanorod growth 
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We study in this paper the atomic mechanisms of nanorod growth and propose the way of diameter 
selection of nanorod. A characteristic radius is demonstrated to be crucial in nanorod growth, which 
increases proportional to one fifth power of the ratio of the interlayer hopping rate of adatoms across 
the monolayer steps to the deposition rate. When the radius of the initial island is larger than 
this characteristic radius, the growth morphology evolves from a taper-like structure to a nanorod 
with radius equal to the characteristic radius after some transient layers. Otherwise the nanorod 
morphology can be maintained during the growth, with stable radius being limited by both the 
radius of the initial island and the three-dimensional Ehrlich-Schwoebel barrier. Therefore different 
growth modes and diameter of nanorod can be selected by tuning the characteristic radius. The 
theoretical predictions are in good agreement with experimental observations of ZnO growth. 

PACS numbers: 61.46.-w, 68.65.Ac, 68.55.-a 



I. INTRODUCTION 

One-dimensional nanostructures have attracted much 
attention since they provide potential applications for 
nanoelectronics and nanophotonics [TH!]. For example, 
zinc oxide nanorods with hexagonal cross-section can be 
applied as whispering gallery resonators, in which the 
coupling between the resonant modes and free excitons 
depends sensitively on the cross-sectional radius [3H5]. 
Growth of nanorods have been extensively reported thus 
far, yet there are few studies concentrating on the under- 
lying atomistic mechanisms of growth, especially on the 
understanding and controlling the cross-sectional radius 
of nanorods from atomic point of view [61-113]. 

It is known that surface kinetics plays an important 
role in determining the morphology and size of nanostruc- 
tures [2HI7]. The deposited adatoms can either diffuse 
within the topmost layer and aggregate to form a new 
layer nucleus, or hop downward across the step edges and 
contribute to the lateral growth of the topmost layer. Un- 
der a certain deposition condition, the kinetics-controlled 
competition between the growth in the normal direction 
to the substrate and the lateral growth is expected to 
determine the growth modes and morphologies [151 116] . 

The surface kinetics can be described by intralayer 
and interlayer hopping rates of adatoms, v = 

exp(— Ed/kT) and v' — v' exp(— E s /kT), respectively, 
where the prefactors vq and u' Q are the attempt rates 
which are approximately of the same value; k is the 
Boltzmann's constant and T the temperature. The in- 
terlayer diffusion barrier (E s ) is normally larger than the 
intralayer one (Ed). The difference of these two values 
is denoted as E es , which is known as Ehrlich-Schwoebel 
barrier (ESB) p~8l[19]. The ratio of v/v' increases with 
ESB as exp(i? es /fcT). The ESB is reported to increase 
with the step height and saturate in several atomic layers, 



with a value usually referred as three-dimensional (3D) 
ESB. The conventional ESB is applicable for a mono- 
layer step and is hereafter termed as two-dimensional 
(2D) ESB [201 [ST]. 

When ESB is small enough to allow sufficient inter- 
layer diffusion, layer-by-layer growth occurs, whereas 
larger ESB leads to multilayer growth [T5]- In the latter 
case, islands initiate from each nucleus, which approach 
nanorods if the cross-sections remain the same along the 
longitudinal direction. The key point to select and con- 
trol the nanorod growth is to understand how the radius 
of the cross-section varies when atomic layers add up un- 
der specific growth conditions. Since there is a large dif- 
ference between the 2D- and 3D- ESBs, it is also crucial 
to identify the specific roles of 2D- and 3D-ESBs in the 
nanorod growth. 

In this paper, we study the influence of deposition rate 
and ESB on the growth process of nanorods. A charac- 
teristic radius has been identified, which increases pro- 
portional to one fifth power of the ratio of the 2D-ESB 
limited hopping rate of adatoms to the deposition rate. 
Both the growth modes and the nanorod diameter can be 
selected by tuning this characteristic radius. We demon- 
strate that when the radius of the initial island is larger 
than the characteristic radius, the growth morphology 
evolves from a taper-like structure to a nanorod with ra- 
dius equal to the characteristic radius. However, if the 
characteristic radius becomes larger than the radius of 
the initial island, by increasing the 2D-ESB limited hop- 
ping rate or decreasing the deposition rate, nanorod mor- 
phology can be maintained during the growth, with a sta- 
ble radius being limited by both the radius of the initial 
island and 3D-ESB. The theoretical predictions of the 
characteristic radius are demonstrated with experimen- 
tal observations of ZnO growth, and good consistency 
has been found. 
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II. NUCLEATION ON TOP OF AN ISLAND 

Let us consider a nanostructure with thickness of n 
atomic layers. For simplicity, the cross section is taken 
as circular. The radius of the i-th layer is denoted as 
R il measured in terms of the surface cell parameter ao- 
The dimensionless area and perimeter are Ai = wR? and 
Li = 2-nRi, respectively. 

Assuming the growth units are deposited in the nor- 
mal direction of the surface, at rate F per surface cell 
of area ciq. The number of adatom 77 per surface cell is 
determined by the diffusion euqation, drj/dt = vV 2 rj + F. 
By solving the diffusion equation, the distribution of the 
number density of adatom can be obtained, 



V 



F 



(R 2 n -r 2 )+ Ve , 



(1) 



where rj e is the dimensionless number density of adatom 
at the edge of A n . Before nucleation occurs on top of 
A n , a deposited adatom on A n has no other choice but 
to hop across the step edges after a survival time of t. 
At steady state, the number of atoms leaving the surface 
per unit time, L n r] e v' , is balanced by that of the atoms 
deposited on the surface per unit time, FA n . It gives the 
number density of adatom on the boundary, 



r) e = {FR n )/{2v'). 



(2) 



The total number of adatom on A n can be obtained by 
integrating, 
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The average survival time of an adatom is thus r = N At, 
where At = l/(FA n ) is the time interval between subse- 
quent deposition events on A n . In addition to At and r, 
another concerned time scale is the traversal time for an 
atom to visit all the sites of A n , r tr = A n /v. 

As A n grows, the probability of nucleation on A n in- 
creases. Once a new nucleus forms on A n , the number 
of atomic layers n increases by one, which leads to the 
growth in the normal direction to the substrate. For the 
simplest case that a dimer is the smallest stable island, 
the nucleation rate on A n can be given by Q = P1P2/ At, 
where p\ = 1 — exp(— t/ At) is the probability that an 
atom is deposited during the presence of another atom 
on the surface, andp 2 = 1 — cx p(— r / T tr) is the encounter 
probability [2"2] . 

It has been reported that with slow deposition the total 
number of adatoms is usually much less than unity, i.e. 
N < 1, which means that r < At [16j[22]. Furthermore, 
in typical island growth with large ESB, v jv' is much 
larger than dimensionless R n , so the first term in Eq. ([!]) 
is much smaller than the second term, and the number 
density of adatom on A n is approximately uniform, r\ ~ 



The total number of adatom becomes N 



which gives 
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(4) 



It indicates that t 3> r tr when vjv' 3> R n , and the 
encounter probability of two adatoms simultaneously 
present on the island pi is approximately unity. The 
nucleation rate therefore can be approximated as 



n = 



F 2 Al F 2 n 2 R 5 r 



L n v' 
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(5) 



Depending on the height of the steps across which the 
adatoms hop to the lower layer, 2D-ESB or 3D-ESB plays 
the role to influence the interlayer atomic diffusion, re- 
spectively. Correspondingly subscripts 2D and 3D will 
be added to v' or £1 in the context to show such a differ- 
ence. 



III. CHARACTERISTIC RADIUS 

For a buried layer, i.e. the atomic layer above which 
a second layer has formed, the condition r 3> r tr means 
that an adatom can always be trapped by the ascending 
steps before getting chance to hop to the lower layer. The 
adatoms on A n before second-layer nucleation, however, 
has no other choice but to hop to A n _i. The lateral 
growth of the topmost layer A n is thus contributed by 
the atoms deposited on area A n _i, 



dAn/dt = FA n _i 



FttR 2 _ 



(6) 



We assume for the moment that ^4„_i is large enough so 
that A n is always smaller than A n _\. 

The probability / that a second layer has nucleated on 
A n , according to df jdt = — /), increases with A n as 
the following, 



/ = 1 - exp(-J), 



where 



r An dt 
1=1 n 2D ^rdA 
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dA 
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(8) 



can be regarded as the average number of nucleus on A„ 
Substituting Eqs. (5) and ^ into Eq.g, 
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where R c is defined as 
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FIG. 1: Schematic of the two different growth modes from the 
initial island (grey colors) with radius of Ri: (a) R\ > R c ; 
(b) Ri < R c - Ro denotes the average radius of the substrate 
occupied by per island, and R n is the radius of the topmost 
layer in the nanorod with n grown atomic layers. 



v' 2D in Eq. ( 10 ) denotes the hopping rate of adatom from 
A n to A n —i, with subscript 2D emphasizing that R c is 
determined by the 2D-ESB across the monolayer step 
edges. Since the probability / goes rapidly from nearly 
zero to nearly unity when 1=1, the condition that 1=1 
can be used as a criterion for the formation of the (n + 
l)-th layer [15]. According to Eq. (pi), the second layer 
nucleus has formed before R n reaches R n -\ if R n -i > 
R c . If R n -i < R c , however, the probability that second- 
nucleus forms on top of A n is still nearly zero even though 
R n reaches R n -i- 

As R c is determined just by the 2D-ESB and the 
growth conditions, such as the deposition rate and the 
growth temperature, it can be regarded as a characteris- 
tic radius of the growth system. In heterogeneous growth, 
it is known that the effect of the foreign substrate on sur- 
face kinetics properties depends strongly on the thickness 
of the grown layers. The interlayer hopping rate u' 2D is 
therefore variable, especially in the first two layers. Ac- 
cordingly, we denote hereafter the characteristic radius 
of the first layer and the other layers as R c q and R c , 
respectively, in order to distinguish their difference. 

It is known that R c q is critical in determination of 
the growth mode, such as layer-by-layer growth or island 
growth in the beginning of heterogenous growth. Here 
we propose that, once the island growth sets in, it is 
the characteristic radius R c that plays a key role during 
the development of a separate island in selecting growth 
mechanisms and the lateral size of the island. 



IV. TWO GROWTH SCENARIOS 

In heterogenous growth, the average radius of the 
foreign substrate occupied by each nucleus is denoted 
as Rq- If Ro is larger than R c o, a second-layer nu- 
cleus forms when the radius of the first layer approaches 
R\ = R c o (Ro/ Rco) 2 ^ 7 < Ro, which means that the sec- 
ond layer nucleus forms before the first layers coalescence 
and thus island growth sets in. We define i?i as the ra- 
dius of the initial island from which the island growth 
starts. 



As discussed in previous section, in heterogenous 
growth the characteristic radius for the second layer 
changes from R c o to R c . In homogeneous growth, or 
in late stage of the heterogenous growth where substrate 
effect is negligible, the characteristic radius can vary by 
changing the deposition rate F or the temperature T, 
according to Eq. (10). In these cases, R c o and R c are de- 



fined as the characteristic radius before and after chang- 
ing the growth conditions respectively, Rq and i?i cor- 
respond to the radius of the topmost two layers at the 
moment that a nucleus forms on Ri after changing the 
growth conditions. 

Two growth scenarios can be identified according to 
the way that the characteristic radius changes. The 
first scenario occurs when R c decreases from R c o, i.e., 
R c < R c o- Since i?i > R c o, it is still larger than R c . 
Thus the third layer forms atop when i?2 approaches 
R c (Ri/ R c ) 2 ^ 7 . If we assume all the buried layers cease 
growing, the radius of each layer is fixed at the moment 
when it is buried by a new layer. Therefore the radius 
of the z-th layer in a nanostructure with n grown atomic 
layers, Ri, is determined by setting I(Ri) = 1, 



Ri 

Rc 



Ri- 



Rc 



Ri 
Rc 



,(i = 2,3,...,n), (11) 



Eq. (11) indicates that Ri decreases rapidly with in- 
creasing i, until it approaches R c . Correspondingly, 
the growth morphology changes from a tapered one to 
a nanorod in several transient layers, as schematically 
shown in Fig. 1(a). 

Obviously this is merely a limiting case in which some 
strong screening effects exist so that the topmost layer 
dominates the deeper layers in capturing the deposit 
atoms [23] ■ The opposite limiting case is that in which 
the growth units are equally deposited on the exposed 
area, thus the buried layers can also grow. In the latter 
case island growth leads to the well-known wedding-cake 
morphology |16| . 

The realistic situation is that between these two cases, 
in which only some finite topmost layers are involved 
in capturing the deposit atoms as a result of a medi- 
ate screening effect. To illustrate the screening effect on 
the island growth, we consider that only the topmost fi- 
nite N g layers keep growing laterally. The quantity 1/N g , 
which is in the range of to 1, can be taken as a measure 
of the screening strength. We have carried out numerical 
calculations of the rate equations with different values 
of N g , the details of which will be reported elsewhere. 
We find that when the number of the atomic layers of 
the island n is smaller than N g , the island grows with 
the well-known wedding-cake shape. When n increases 
larger than N g , the radii Ri for i < N g grow gradually to 



Ro, while for N g < i < n 



N g Ri approach their stable 
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values after sufficient growth, 



Rr. 



(N g <i<n-N g ). 



(12) 



For 



X is a constant determined by (Ri/R c ) and Ng 
N g = 1, X = (i?i/i? c ) 35 , consistent with Eq. (fl). The 
radius Ri decreases with i, until it approaches R c af- 
ter some transient layers and then remains this value. 
The number of the transient layers is proportional to Ng , 
with coefficient a of the order of magnitude of 10. We 
therefore show that, under a certain screening strength, 
wedding-cake morphologies (n < N g ), tapered morpholo- 
gies (N g < n < aNg) and nanorods (n > aNT) can be 
observed successively during the island growth. 

It is thus clear that if the radius of the initial island 
R\ > i? c , the structure finally approaches a nanorod 
with a tapered base beneath. The screening strength only 
influences the number of the transient layers of conver- 
gence, i.e., the atomic layers of the tapered base. Since 
the nucleation takes place on the topmost monolayer, 
we refer the nanorod converged from Ri > R c as the 
2D-ESB-limited nanorod, with radius R2D equals to the 
characteristic radius of the system, R 2 d = R c - We term 
this growth mode as the 2D-ESB-limited one. 

The second scenario occurs when R c increases from R c q 
to a value much larger than R± . In this case the average 
number of nucleus on top of A2 when A2 covers A\ is 
(Ri/R c ) 5 ~ 0. Therefore the topmost two layers can 
bunch to a bilayer, which then grows laterally from Ri 
to i?2 till a new nucleus finally forms atop. The process 
repeats and the growth morphology remains rod-like as 
shown in Fig. 1(b). The whole nanorod grows laterally 
as the number of the atomic layers n increases, fed by the 
deposited adatom on the top of the nanorod. Therefore, 



— ^ = —FA n = —F-kR 
at n n 



(13) 



The average number of nucleus can be obtained by inte- 
grating Eq. (HI, 
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where R, 
Eq, 



is the same characteristic radius defined in 
(10). Note H,3D denotes the nucleation rate on A n 
after A n approaches A n -i, and v' 3D represents the in- 
terlayer hopping rate across the multilayer step edges. 
As discussed above, when R c is much larger than R n -i, 
the first term in Eq. ( 14 1 is nearly zero. The radius of 
nanorod with n atomic layers R n can be therefore deter- 
mined by setting / 



1 in Eq. (14), 
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FIG. 2: The radius of the n-th layer R n in an island with n 
grown atomic layers, for the cases of i?i > R c with N g — 1 
(above the dotted line), and Ri < R c (below the dotted line), 
as a funct ion of n. The values are calculated according to Eqs. 
(Ill and ( 16 1. The dotted line guides the characteristic radius 
R c - a = v 3D /v' 2D is the ratio of the hopping rate across the 
multilayer step to the one across the monolayer step. 



where a 



, 3 M D = cxp[-(£ e 3 f - E™)/kT] and < 



a < 1. The radii of the i — th layer Ri can be obtained 



by iterating Eq. (15), 



Ri Rn 
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The nanorod radius increases until the difference of 
R n and R n -i is smaller than one lattice parameter for 
sufficient large n, when it approaches the stable radius 
of the 3D-ESB-limited nanorod R^d- According to Eq. 



( 15 1, this happens when 



(17) 



Substituting Eq. (17) into Eq. (16) and replacing the 



summation with logarithm for large enough n, the stable 
radius of the 3D-ESB-limitcd nanorod can be obtained 



as, 



R3D 



R\ + -aRl( n 



1 + ln 



aRl . 

lri 3D 



1/5 



(18) 



wheren 7 is the Euler's constant. Since aR^, is propor- 
tional to v' 3D /F, it is evident that the nanorod growth in 
this case is determined by the radius of the initial island 
Ri, the deposition rate and the 3D-ESB. Larger 3D-ESB 
( smaller a) facilitates the convergence (i.e. smaller n s ) 
and leads to smaller R 3 d, which is consistent with a re- 
cent Monte-Carlo simulation on copper nanorod growth 



It is worthy to emphasize that Eq. (15) is valid only 



when the first term in Eq. ( 14 ) is negligible, and it is 



(15) 

physically invalid to extrapolate from Eq. (|18|) to a R 3 d 
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FIG. 3: (a) Schematic of the temperature change in exper- 
iments, where the three dashed lines correspond to the mo- 
ments the growth is cut off by large flux of nitrogen. The ZnO 
products obtained at the three cut off moments are shown in 
(b), (c) and (d), respectively. 



larger than R c . Actually, the 2D-ESB-limited growth 
sets in once R n is increased to R c . Therefore the real 
radius R^d may never exceed R c . 

For comparison, we show in Fig. 2 the radius R n in a 
nanostructure with n grown atomic layers for the two sce- 
narios, according to Eqs. (Ill and (16 1. It is clear that 



when the radius of the initial island R\ is larger than 
the characteristic radius R c , R„ decreases rapidly until 
it approaches R c . Consequently, the growth morphology 
evolves from a taper-like structure to a nanorod with uni- 
form radius of R c . In this scenario, the converged radius 
is limited by 2D-ESB. When R\ < R c , Rn corresponds to 
the nanorod radius. It increases relatively slowly, with a 
stable radius smaller than R c , determined by the 3D-ESB 
and the radius of the initial island R\ . 



V. EXPERIMENT VERIFICATIONS 

Experimentally we take zinc oxide (ZnO) vapor growth 
as an example to verify the selection of the nanorod ra- 
dius via varying R c by changing the growth tempera- 
ture. For this purpose, unlike conventional ZnO nanorod 
growth systems where catalysts are usually introduced, 
here we establish a physical growth system without using 
additive chemicals. 

The nanorods of ZnO were synthesized catalyst-free in 
a horizontal tube furnace with programmable tempera- 
ture control. The pure zinc powder (99.9% Alfa Aesar) 
and polished Si(100) substrate were arranged in the same 
quartz boat and 1.0 cm apart. The growth was carried 
out with flux of nitrogen controlled as 300 standard cu- 
bic centimeter per minute (seem) and oxygen gas flux as 
5 seem. The temperature in the central section of the 



furnace was homogeneous, where the quartz boat was 
placed. In each run of the experiment, the temperature 
was changed as shown in Fig. 3(a) to control the depo- 
sition rate. The growth of nanorods was terminated at 
different time by sudden increasing the nitrogen flux and 
cutting off the oxygen supply, as illustrated by the three 
dashed lines in Fig. 3(a). In this way, we expect to pre- 
serve the growth morphology at the moment that the 
growth has been terminated. 

In our experiments, the temperature is varied while 
the flux of N2 and O2 are kept as constants for ZnO 
nanorod growth. The evaporated zinc atoms react with 
oxygen molecules, and the partial pressure of the produc- 
tion ZnO is proportional to that of zinc, PznO cx Pzn/K p , 
where Kp is the reaction constant. Both P% n and K p ex- 
ponentially depend on temperature, 



P Zn oc exp(- — ), 



B 



if P cxcxp(~-^), (19) 



where B Zn = 0.58 eV (6776K) [H] and B K = 0.21 eV 
(2474 K) [26]. The partial pressure of ZnO can be there- 
fore written as, 



Pzno = Po eM-B/kT) 



(20) 



where B — B% n — B^ = 0.37eV, and Po is a constant de- 
termined by the other growth conditions except of tem- 
perature. The temperature-dependent deposition rate 
per lattice site can be written as F — eto-PznO / 'v lc lixmkT '. 
According to Eq. (10), the characteristic radius R c is 



R c = sb(kT) 1/10 exp(-AE/kT), 



(21) 



where b = [7u V2Tvm/(P Q al)] 1 , AE = (E s -B)/5, and 
s is the geometrical factor associated with the different 
cross-sectional shapes of nanorod. Equation (21 1 shows 



obviously that R c can be tuned by changing temperature. 
If AE is positive, R c decreases with decreasing temper- 
ature, then the first (2D-ESB limited) growth scenario 
is realized, and the radius of nanorod corresponds to the 
characteristic radius R c . Otherwise if AE is negative, R c 
increases with decreasing temperature, and so does the 
radius of nanorods. 

Now let us check the variation of morphologies of ZnO 
nanorods with step-decreased temperature as shown in 
Fig. 3(a). The morphologies of the nanorods were charac- 
terized by a field emission scanning electron microscopy 
(LEO 1530VP). Figure 3(b) illustrates ZnO structures 
grown at constant temperature of 600° C. There is no ev- 
ident variation of the cross-sectional diameter along the 
longitudinal direction of the nanorods. When the tem- 
perature was decreased from 600 to 550 °C, a second 
segment of nanorods appear, the cross-section area of 
which shrink to smaller values, as shown in Fig. 3(c). The 
morphologies obtained after two temperature drops are 
shown in Fig. 3(d), where two evident changes of cross- 
sectional diameter can be identified on the nanorods as 
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FIG. 4: The diameter of the second segment nanorod as a 
function of the corresponding temperature. The dashed curve 
gives the theoretical fitting results according to Eq. (21 1. 



dius is 2D-ESB-limited and approaches the characteristic 
value. Otherwise the nanorod is 3D-ESB limited with a 
stable radius smaller than the characteristic radius. We 
suggest that our results is helpful to select the desired 
growth modes and control the diameter of nanorods and 
nanowires. 

Although experimental studies on growth of ZnO 
nanorod with hexagonal cross-section has been reported 
before, to the best of our knowledge, a quantitative study 
considering the kinetics in the interfacial growth remains 
rare. Moreover, the theoretical model proposed here in 
fact is a generic one that is not limited to ZnO nanorod 
growth only. Experimentally, if one can precisely tune 
the characteristic radius R c for a specific growth system, 
or choose a desired initial radius by using suitable seeds, 
either kinds of growth modes can be selected in order to 
obtain different morphologies and nanorod size. By this 
means, microscopic informations of 2D-ESB and 3D-ESB 
can also be inferred from the experiments. 
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structures as shown in Fig. 3(c), as a function of the 
corresponding temperature. The dashed curve gives the 
theoretical fitting results according to Eq. (21). It shows 



that the theoretical model is in good consistency with 
the experimental data. The fit value of AE is 0.59 eV, 
which leads to E s of about 3.3 eV. We should point out 
that this value is only a rough estimate of the effective 
barrier against the adatom diffusing across a monolayer 
step edge in zinc oxide system. 



VI. SUMMARY 

We have demonstrated that two nanorod growth 
modes can be realized, depending on a characteristic ra- 
dius which is determined by the 2D-ESB and the de- 
position rate. When the radius of the initial island is 
larger than this characteristic radius, the nanorod ra- 



This work was supported by NSF of China (10974079 
and 10874068) and Jiangsu Province (BK2008012), 
MOST of China (2004CB619005 and 2006CB921804). Z. 
Zhang acknowledges partial support by USDOE (grant 
No. DE-FG02-05ER46209, the Division of Materials Sci- 
ences and Engineering, Office of Basic Energy Sciences) , 
and USNSF grant No. DMR-0906025. 



* Electronic address: djshu@nju.edu.cn 
' Electronic address: muwang@nju.edu.cn 
[1] S. V. N. T. Kuchibhatla, A. S. Karakoti, D. Bera, and 

S. Seal, Prog. Mater Sci. 52, 699 (2007). 
[2] M. H. Huang, S. Mao, H. Feick, H. Q. Yan, Y. Y. Wu, 

H. Kind, E. Weber, R. Russo, and P. D. Yang, Science 

292, 1897 (2001). 
[3] T. Nobis, E. M. Kaidashev, A. Rahm, M. Lorenz, and 

M. Grundmann, Phys. Rev. Lett. 93, 103903 (2004). 
[4] L. K. van Vugt, S. Riihle, P. Ravindran, H. C. Gerritsen, 

L. Kuipers, and D. Vanmaekelbergh, Phys. Rev. Lett. 

97, 147401 (2006). 
[5] L. Sun, Z. Chen, Q. Ren, K. Yu, L. Bai, W. Zhou, 

H. Xiong, Z. Q. Zhu, and X. Shen, Phys. Rev. Lett. 100, 

156403 (2008). 

[6] X. Wang, C. J. Summers, and Z. L. Wang, Nano Lett. 4, 
423 (2004). 

[7] J. Y. Lao, J. G. Wen, and Z. F. Ren, Nano Lett. 2, 1287 
(2002). 

[8] B. H. Kong, D. C. Kim, and H. K. Cho, Physica B 376, 
726 (2006). 

[9] K. Kitamura, T. Yatsui, M. Ohtsu, and G. C. Yi, Nan- 
otechnology 19, 175305 (2008). 
[10] J. Xiao, Y. Wu, X. Bai, W. Zhang, and L. G. Yu, J. Phys. 
D: Appl. Phys. 41, 135409 (2008). 



7 



[11] R. C. Wang, C. P. Liu, J. L. Huang, S.-J. Chen, Y.-K. 
Tseng, and S.-C. Kung, Appl. Phys. Lett. 87, 013110 

(2005) . 

[12] A. Umar and Y. B. Hahn, Appl. Phys. Lett. 88, 173120 

(2006) . 

[13] Z. W. Liu and C. K. Ong, Mater. Lett. 61, 3329 (2007). 
[14] Z. Y. Zhang and M. G. Lagally, Science 276, 377 (1997). 
[15] J. Tersoff, A. W. Denier van der Gon, and R. M. Tromp, 

Phys. Rev. Lett. 72, 266 (1994). 
[16] J. Krug, P. Politi, and T. Michely, Phys. Rev. B 61, 

14037 (2000). 

[17] P. Politi, G. Grenet, A. Marty, A. Ponchet, and J. Villain, 

Phys. Rep. 324, 271 (2000). 
[18] G. Ehrlich and F. G. Hudda, J. Chem. Phys. 44, 1039 

(1966). 

[19] R. L. Schwoebel and E. J. Shipsey, J. Appl. Phys. 37, 



3682 (1966). 

[20] S. J. Liu, H. Huang, and C. H. Woo, Appl. Phys. Lett. 

80, 3295 (2002). 
[21] M. G. Lagally and Z. Y. Zhang, Nature 417, 907 (2002). 
[22] J. Rottler and P. Maass, Phys. Rev. Lett. 83, 3490 

(1999). 

[23] S. van Dijken, L. C. Jorritsma, and B. Poelsema, Phys. 

Rev. Lett. 82, 4038 (1999). 
[24] L. G. Zhou and H. Huang, Phys. Rev. Lett. 101, 266102 

(2008). 

[25] C. B. Alcock, in CRC Handbook of Chemistry and 
Physics, edited by D. R. Lide (CRC Press, Boca Raton, 
FL, 2004), 84th ed. 

[26] Z. L. Wang, J. Phys.: Condens. Matter 16, R829 (2004). 



